perm filename SBOUND.CH[UHF,DEK] blob
sn#841772 filedate 1987-06-19 generic text, type T, neo UTF8
% change file to compute the error bounds and actual errors
@x
if k<62 then
begin err←err*alpha[k+2];
for l←start[k] to start[k+1]-1 do
begin u←i+del_i[l]; v←j+del_j[l];
buffer[u,v]←buffer[u,v]+err;
@y
if buffer[i,j]<min_d then
begin min_d←buffer[i,j]; write_ln(tty,min_d:10:5);
end
else if buffer[i,j]>max_d then
begin max_d←buffer[i,j]; write_ln(tty,max_d:10:5);
end;
if k≥62 then
begin tot_n←tot_n+1; tot_err←tot_err+abs(err);
if abs(err)>1.0 then tot_big←tot_big+1;
end
else begin err←err*alpha[k+2];
for l←start[k] to start[k+1]-1 do
begin u←i+del_i[l]; v←j+del_j[l];
buffer[u,v]←buffer[u,v]+err;
if buffer[u,v]<min_d then
begin min_d←buffer[u,v]; write_ln(tty,min_d:10:5);
end
else if buffer[u,v]>max_d then
begin max_d←buffer[u,v]; write_ln(tty,max_d:10:5);
end;
@z
@x
@!alpha:array[0..63] of real; {multipliers and thresholds}
@y
@!alpha:array[0..63] of real; {multipliers and thresholds}
@!tot_err:real; {accumulated absolute value of errors}
@!tot_n,@!tot_big:integer; {number of barons and number of times |err>0.5|}
@!min_d,@!max_d:real; {min and max density so far}
@ @<Set init...@>= tot_err←0; tot_n←-126; tot_big←0; min_d←0.0; max_d←0.0;
@z
@x
write_ln('\endddtone');
@y
write_ln('\endddtone');
write_ln(tty,tot_n:1,' barons; ',tot_big:1,' bigs; ',tot_err/tot_n:10:5);
@z